Equilibration of an isolated quantum system and thermalization of its subsystem 
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An isolated quantum chaotic system in a pure state in general can equilibrate and serve as a heat 
bath. We show that once the equilibrium is reached, any of its subsystems that is much smaller 
than the isolated system is thermalized so that the subsystem is governed by the Gibbs distribution. 
Within this theoretical framework, the celebrated superposition principle of quantum mechanics 
should lead to a prediction of a thermalized subsystem with multiple temperatures if the isolated 
system is in a superposition state of energy eigenstates of multiple distinct energy scales. This 
multiple temperature state is at equilibrium, completely different from a non-equilibrium state that 
has multiple temperatures at different parts. Possible experimental schemes to verify this prediction 
are discussed. 



Standard textbook always starts its discussion of quan- 
tum statistical mechanics with the microcanonical en- 
semble where every energy eigenstate in a narrow en- 
ergy range is equally possible. Mathematically, this 
ensemble is described by the density matrix p mc = 
\a n \ 2 \E n )(E n \ , where E n lies in a narrow energy 
range [E, E + AE] and |a„|'s are equal to each other. Mi- 
crocanonical ensemble usually refers to an isolated sys- 
tem and should be described by a wave function. To 
derive the density matrix p mc from a wave function, 
two postulates-equal a priori probability and random 
phases-have to be introduced. In standard textbooks, 
these two postulates are either explicitly stated pQ or 
implicitly assumed [2]. One way to understand these two 
postulates is to imagine that the microcanonical ensem- 
ble is not truly isolated but weakly interacting with an 
outside system pQ. This seems reasonable as a system 
that we study always interacts with an external world. 
However, one is always allowed to expand the system to 
include more systems and eventually have a truly isolated 
system. As a result, it is legitimate to ask whether we can 
derive the main results of quantum statistical mechanics 
starting with a truly isolated system[3]. An equivalent 
question is whether we can derive thermodynamics from 
the dynamic evolution of a wave function [3] . 

There have been increased efforts to address the issue 
of equilibration of an isolated quantum system [^I30j. 
partly due to that an isolated quantum system can now 
be achieved and sustained experimentally for a reason- 
able long time in a highly excited state with ultra-cold 
atoms [3Tj . Much progress has been achieved through 
these efforts. For example, Rigol et al. demonstrated 
numerically that the momentum of an interacting bo- 
son system in a highly excited eigenstate satisfies the 
microcanonical prediction, giving strong support to the 
eigenstate thermalization hypothesis (ETH) (5TU0|. It 
was proved analytically that an isolated quantum sys- 
tem satisfying certain loose conditions can indeed equili- 
brate [HI H UllSI- 
In this work we focus on quantum chaotic systems [32j- 
IM] to address the issue of equilibration of an isolated sys- 



tem. We show numerically with the Henon-Heiles system 
that a general quantum state in a chaotic system does 
equilibrate after a long time evolution. This result is not 
only consistent with an inequality for equilibration of an 
isolated system obtained previously [HHHj, but also pro- 
vides an intuitive picture on how this equilibration takes 
place. We then prove analytically that a subsystem of 
an isolated quantum system at equilibrium is thermal- 
ized such that it is described by the Gibbs distribution. 
Here we distinguish between equilibration and thermal- 
ization: a system equilibrates if its overall features no 
longer change with time; a system is thermalized if its 
state is described by the Gibbs distribution. An equili- 
brated system is not necessarily thermalized while a ther- 
malized system must be at equilibrium. 

A natural and surprising outcome of this theoretical 
framework is that a subsystem can thermalize with mul- 
tiple distinct temperatures. This can happen when the 
isolated system is in a superposition of energy eigen- 
states that concentrate around different energy scales. 
Although to the best of our knowledge it has never been 
observed in nature and in experiment, a thermalized sys- 
tem with multiple temperatures appears unavoidable as 
there is no a priori reason that an isolated large quan- 
tum system must be in a state which is composed only of 
energy eigenstates from a narrow energy range. We em- 
phasize that this multiple-temperature state is at equilib- 
rium so that it is completely different from a state that 
is out of equilibrium and has different temperatures at 
different parts of the system or for its different degrees of 
freedom. Possible experimental schemes to confirm our 
predication are discussed. 

Equilibration of an isolated quantum chaotic system. — 
An important result has been achieved for the equili- 
bration of an isolated quantum system. By generalizing 
Reimann's result [TS] , Short et al. [TTJ [T2J proved an in- 
equality regarding the fluctuation of an observable A in 
the course of time evolution 

o\ = (\tr{A \m) W)\} tr(A Pco )| 2 ) t < ^ , (1) 
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FIG. 1: Time evolution of a wave-packet and the long-time average in the Henon-Heiles system. The first row is the density in 
the real space and the second row is the density in the momentum space. The long-time averages of these densities are shown 
at the rightmost panels. The average is taken over 1200 states equally separated in the time interval of [0.2012,0.2255]. The 
unit for the real space is r c and the unit for the momentum space is po- The red lines in the first row are energy contours 
of the Henon-Heiles potential at V(x,y)/V c = 1/2, 1, 2. The red lines in the second row are the maximal classically allowed 
momentum for the initial energy of the Gaussian wave packet. The color bars are given on the right side. 



where p{t) = \^{t)) (^{t)\ with = E k c k \E k ) be- 

ing the wave function of the isolated system and p^ — 
12k \ Ck \ 2 \Ek) (Ek\- The subscript t in () t represents an 
averaging over a time period much longer than the char- 
acteristic time scale of the system. ||A|| 2 is the maximal 
value of A regarding all the states in the Hilbert space. 
The effective dimension d„a = „ , = * . indi- 

cates how widely the state |^) is spread over the energy 
eigenstates. For the micocanonical ensemble density ma- 
trix p mc , one can easily check that d e Q = N with N being 
the number of energy eigenstates occupied. This inequal- 
ity holds for a large class of quantum systems, including 
quantum chaotic systems. 

If a quantum system is in a state of high energy, then 
d e Q should be large since the density of states usually 
increases tremendously with energy. This means that 
the right hand side of Eq.([T]) is small. For this case, 
this inequality tells us two things: (i) An isolated quan- 
tum chaotic system in a high-energy state will eventually 
equilibrate so that a bounded observable will fluctuate in 
small amplitude around its averaged value for almost all 
times. By almost all times, we mean that the fluctua- 
tion does get very large occasionally. This is similar to 
the familiar fact in classical statistical mechanics: even 
when the equilibrium is reached, large fluctuations can 
still occur rarely, (it) Although the isolated quantum 
chaotic system is always described by a wave function, 
the expectation value of an observable A at almost any 
moment can be computed with density matrix p^,, that 
is, tr(|?^) (ip\ A) tr(p 00 A). Note that density matrix 
Poo is different from the microcanonical density matrix 



p cm as |cfc| 2 's are determined by the initial condition and 
they are not necessarily equal to each other. 

The inequality in Eq.Q dictates that an isolated 
quantum chaotic system in high energy state can equi- 
librate; however, it does not show how this equilib- 
rium is reached in a dynamical process. We now pro- 
ceed to illustrate numerically this dynamical equilibra- 
tion process with the Henon-Heiles system, a well-known 
chaotic system 35J. The Hamiltonian of the system is 
H = p 2 /2m + f (x 2 + y 2 ) + X(x 2 y ~ The Henon- 
Heiles system has three saddle points located at a dis- 
tance r c = 4 from the origin. These three points are the 
corners of the energy triangular contour with potential 
V c = |r^2 as shown in Figjl|al-a5). The momentum cor- 
responding to the saddle point energy is po = \j2mV c as 
indicated by the circle in Fig{l|bl-b5). In our numerical 
simulation we set m = |, H = 1, and A/X = 1/3. 

The initial condition is a highly localized Gaussian 
wave packet as shown in Figjl] (al,bl) so that the sys- 
tem energy is high. This wave packet is centered at 
fi = (0.3,0)r c and ft = (cos 10°, sin 10°) ^7/10p in the 
real and momentum spaces, respectively. A classical par- 
ticle with ?i and Pi has energy 0.9691^. The motion of 
this classical particle is fully chaotic but confined in the 
triangular region and the circled region in the real and 
momentum space, respectively. 

We numerically solve the time dependent Schrodinger 
equation and the dynamical evolution of the wave packet 
is illustrated in Fig[TJal-a4,bl-b4). As the wave packet 
evolves, it begins to spread out and distort in shape. 
Eventually, it reaches an equilibrium state, where the 
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wave packet spreads out all over the triangular region 
in the real space and the circled region in the momen- 
tum space. This overall feature will no longer change 
even though the details of the wave packet still change in 
the following dynamical evolution. This is similar to the 
equilibrium in classical systems, where the microscopic 
state still changes while the overall macroscopic quan- 
tities such as pressure and temperature do not change. 
This kind of dynamical equilibration process has been 
demonstrated for chaotic billiard systems [HJ [35] . 

For comparison, we have calculated n oc (r) = (r\poo\r) 
and n oc (p) = (P\poo\p) by long-time averaging, i.e. the 
equilibrium state obtained by Short and Rcimann |llLll2l 
[T9] , The results are shown in Fig(lja5,b5). It is clear that 
^oo (f) an d fioo (p) are very similar to the wave packet at 
t = 0.2126 with the same overall feature that the wave 
function spread all over both the triangular spatial region 
and the circled momentum region. The difference is that 
the wave packet at t = 0.2126 has larger fluctuations, 
whose distribution obeys an exponential law |21j while 
the fluctuations in both n^^r) and n^p) are Gaussian. 

To quantitatively illustrate the dynamical relaxation 
process, we follow the idea of Kolmogorov and Sina [57] 
and define an entropy as S$ — j dt; ln( " jl) ) > £ — 
r, p. The time evolution of is shown in Fig|2Tal,a2). 
In the figure, we see clearly the entropies quickly satu- 
rate and reach the maximum values, indicating that an 
equilibrium state is reached. We also notice that the re- 
laxation times in both the real and momentum spaces 
are about the same. 

The equilibrium state reached is consistent with 
the inequality for equilibration obtained by Short and 
Reimann[TTJ [T^l HH]. To check the inequality numer- 
ically, one needs to compute energy eigenstates of the 
system. As it is difficult to compute the eigenstates for 
the Henon-Heiles system, we have turned to the ripple 
billiard system studied in Ref . [STJ [35] , where the inequal- 
ity is verified. We will present the detailed computation 
elsewhere and only mention that d e e ~ 300 for a similar 
Gaussian wave packet in the ripple billiard system. 

We have further explored the link between the quan- 
tum dynamical relaxation discussed here and its classi- 
cal counterpart. For this purpose, we compare density 
matrix with the classical microcanonical ensemble 
p c = ^5(H(p,r) — E) in both the real and momentum 
spaces. For better comparison, we choose the following 
marginal distribution: for the real space, we integrate 
out the y dimension to obtain the density distribution 
P(x)] for the momentum space, we integrate out the 
angle variable to find the momentum distribution f(p). 
The results are shown in Fig{2]jbl,b2), where we see the 
quantum and classical distribution are consistent, except 
some quantum tunneling effect indicated by the non-zero 
value in the classically forbidden region. This link in- 
dicates that the quantum dynamical relaxation process 
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FIG. 2: (a) Time evolution of the Kolmogorov-Sina entropies 
S$ for the Henon-Heiles system for both spatial and momen- 
tum space, (b) Comparison between equilibrium quantum 
and classical distributions in the Henon-Heiles system. P(x) 
is for the real space and f(p) for momentum space. Unit in 
the real space is r c and unit for the momentum space is po. 



can be understood with the help of the classical picture. 
Intuitively, one can think of the wave-packet as a col- 
lection of classical particles, the equilibration of a wave 
packet can then be likened to the relaxation process of 
the classical particles in the phase space. This kind of 
quantum-classical correspondence in dynamics has been 
studied before [3U[ |30~]. However, the previous focus seems 
to be on the motion of the wave packet center; our fo- 
cus is on the relaxation process and the final equilibrium 
state. Finally, we note that the quantum-classical cor- 
respondence shown in Figj2] is not contained in the in- 
equality Eq.Q. There is much more to be explored for 
the quantum dynamical relaxation in the future. 

Thermalization of subsystems — We have shown that a 
truly isolated quantum chaotic system does equilibrate. 
So should any of its subsystems that are small compared 
to the isolated system but large enough to be consid- 
ered as macroscopic objects. We shall show that due 
to the coupling to the rest of the system, these equili- 
brated subsystems are also thermalized so that they are 
described by the Gibbs distribution. The derivation of 
Gibbs distribution for a subsystem has been considered 
before with the assumptoin that the isolated system is in 
a pure state composed of energy eigenstates from a small 
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energy interval [241 128] . We show that this assumption is 
not necessary and when the pure state is composed of en- 
ergy eigenstates of different energy scales, the subsystem 
is thermalized with multiple temperatures. 

We decompose the isolated quantum chaotic system 
into two parts, system S and thermal bath B and 
write the Hamiltonian as H s+B = H s + H B + AH, 
where AH is the weak interaction between system H s 
and thermal bath H B . Suppose that the compos- 
ite system is described by a wave function |$ ,S + B ) = 
J2 c k\E k +B ), where |£ , f +s )'s are the energy eigen- 
states of the composite system. By tracing out ther- 
mal bath B, we obtain the density operator for system 
S, p s = tr (s) \<$> S+B ) ($ S+B |. The system will even- 
tually equilibrate; as an observable, p s will be close to 
its long time average, i.e. p s s» (p s ) t = tr^ B 'p^ B = 
tr^E\c k \ 2 \E S k +B )(E S k +B \. 

We expand the eigenstate \E k +B ) as follows 

l^ +B > - l*f ) \Ef) , (2) 

where \Ef) and \E B ) are energy eigenstates of system 
S and thermal bath B, respectively. The prime above 
indicates that the summation is only over eigen-energies 
satisfying 

E% +B = E? + Ef + AEij . (3) 

where AEij is the interaction energy that is usually very 
small compared to Ef and E B when long-range inter- 
action is negligible, e.g., gravity, in the system. Two 
remarks are warranted, (i) The approximation made in 
Eq.([2]) is justified. The equality holds when there is no 
coupling AH = 0. We expect it hold when the weak 
interaction AH is turned on. (ii) The weak interaction 
AH can drive the system to a state with af-'s randomly 
uniformly distributed on the sphere ^ . |a^| 2 = 1. This 
random distribution is similar to the idea of "typical- 
ity" [22l E3]. The connection between interaction and 
randomness is widely acknowledged since the details of 
the interaction is irrelevant to the statistical proper- 
ties [U H2]- As a result, the average value of |a*-| 2 
is dS+b ^ eS+B) , where D S+B (E^ +B ) is the degeneracy 

brought by the combination of states. We emphasize that 
this degeneracy is different from the intrinsic degeneracy 
of the system and it is due to the existence of AJ3y in 
Eq.@ . 

With the approximation made in Eq.Q, we now pro- 
ceed with our derivation, 

P S = Y.\ c ^Y.^n\Et +B ){E S k +B \E B ) 

k m 

k ii' m 

k i m 



+Ei c fci 2 E(E a ™ a "™)i^fx^i ■ (4) 

k i^V m 

The central limit theorem gives the results of the 
summation £ m |4J 2 * D B (E B n ) / D s + B (E s k +B ), 
where D B (E B l ) is the degeneracy of the ther- 
mal bath. There are more terms in the sec- 
ond summation. However, again due to the 
central limit theorem, the summation £ m a\ m a\* m 
is of the order 0{^/D B (E B )/D s + B (E k s+B )} ~ 

0{l/^J D S + B (E^ +B )}, where we have used that S is 

much smaller than B so that D B (E B ) ~ D S+B (E^ +B ). 
As a result, the last term i n Eq. Q has the order of magni- 
tude at 0{£\ D s (E l ) 2 /^D s + B (E^ +B )}, which is prac- 
tically zero for the isolated system is much larger than 
system S. So omitting the last term, we have from Eq. Q 

P S = EN 2 E osTnfjL i^fl- (5) 

k i \ k ' 

With the standard argument for the Gibbs distribu- 
tion [2], we arrive finally at 

p s = E ki 2 {E ex p(-«) \ E f) ( E "\} < ( 6 ) 

k i 

where (3 k = j^r = d\nD s+B (E^ +B ) / 8E s k +B de- 
fines the temperature of the total system for eigenstate 
In this way, we have proved that a subsystem 
of an isolated equilibrated system is thermalized. 

Thermalized state with multiple temperatures — Here 
we examine Eq. Q for two typical cases: (i) The co- 
efficients |cfc| 2 of the composite system have a single 
sharp peak distribution around energy E p . This case 
is considered by others [24] in different formalisms. For 
this case, the density matrix p s in Eq.([6]) is reduced to 
p s = J2i ex P(-/3p^) iV'f) (ipf\- This is exactly the typ- 
ical Gibbs distribution discussed in all textbook on sta- 
tistical mechanics, (ii) The coefficients \ck\ 2 have two 
well-separated sharp peaks around two energies E v \ and 
E p 2- In other words, the composite system (or the heat 
bath) is in a superposition of numerous eigenstates cen- 
tered around two very different energy scales. In this 
case, the thermalized system has now two temperatures, 
pi for E p i and @2 for E p %, with an unconventional density 
matrix 

P s = E(kl 2 e-^ + N 2 e-^) |J5f) (E?\ , (7) 

i 

where |ai| 2 and | a.2 1 2 are the weight of the two distri- 
bution peaks. This existence of a thermalized system 
with more than one temperatures appears to be an un- 
avoidable result from the superposition principle. When 
the heat bath is quantum, there is no a priori reason that 
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FIG. 3: Momentum distributions of a thermalized ideal gas 
with two temperatures, (a) The two temperatures are of the 
same scale; (b) the two temperatures are of different scales, 
and a pronounced double-peak distribution is seen. Red solid 
line is for the two-temperature distribution; blue dashed line 
for the lower temperature; green dashed line for the larger 
temperature. 



the heat bath must be in the quantum state which is nar- 
rowly peaked at one energy. When the heat bath is in a 
superposition of states with well-separated energy scales, 
its subsystem has to thermalize with multiple tempera- 
tures. 

Here are two caveats, (i) Such a state is an equilibrated 
state with multiple temperatures and it is completely dif- 
ferent from the usual non-equilibrium state which has 
different temperatures for different parts, e.g., in the 
double- well experiment |43j . (ii) The isolated composite 
system is in a superposition of numerous eigenstates cen- 
tered around two very different energies. This is very 
different from a state which is a superposition of only 
two different eigenstates. For the latter case, the system 
does not equilibrate as the right hand side of Eq.Q is 
large with d e s « 2; consequently, its subsystem does not 
thermalize. 

As an example, we have computed the momentum dis- 
tribution f(P) of a three dimensional ideal gas that is 
thermalized with two temperatures. For computation, 
we choose |ai| 2 = 0.3 and | (3,2 1 2 = 0.7 in Eq.Q and have 
both temperatures much higher than the quantum de- 
generate temperature. Two typical results are shown 
in Fig|3] When the two temperatures are of the same 
scale, the two-temperature momentum distribution is not 
much different from a single temperature distribution 
(see Fig{3]ja)); when the two temperatures are of differ- 
ent scale, the two-temperature state has a pronounced 
double-peak distribution (see Fig[3jb)), revealing the ef- 
fect of two energy scales in the quantum heat bath. 

For all the systems that we have encountered in nature 
or studied in experiment, they are in contact with clas- 
sical heat bath. This may explain why such a thermal- 
ized state with multiple temperatures has never been ob- 



served. However, with the advance of technology, we can 
now create large quantum systems such as superfluid he- 
lium and Bose- Einstein condensates (BECs), which may 
serve as quantum heat bath. In addition, it has been 
reported that the quantum and classical fluctuations can 
be separated out with the spin echo technique in a nu- 
clear heat bath [HI H5] . We have also the ability to put 
these systems in highly excited states while keeping their 
coherence [3T] . These technical achievements should pave 
the way for finally observing such an exotic thermalized 
state. One possible experimental scheme is to use the 
two-specie BECs with one specie has much larger pop- 
ulation than the other. One may be able to drive the 
system to have two energy peaks with the double-well 
potential and then observe the double-peak momen- 
tum distribution shown in Figj3]jb) for the smaller specie. 

Note that Fine et al. have also found it not neces- 
sary to assume that a microcanonical ensemble is com- 
posed of energy states from a narrow energy range |46F 
148] - They replaced the standard microcanonical ensem- 
ble with "quantum microcanonical" ensemble, which is 
assumed to consist of all superpositions for a given en- 
ergy expectation. This is fundamentally different from 
our approach, where no assumption for an ensemble is 
needed. However, physically, the quantum microcanoni- 
cal" ensemble considered by Fine et al. may be equivalent 
to a special case in our theory, where the pure state of 
the isolated system is made of equally-weighted energy 
states from a wide energy range. 

Summary — We have shown that a truly isolated quan- 
tum chaotic system can equilibrate and its subsystem can 
thermalize to the Gibbs distribution. Moreover, within 
such a theoretical framework we predict that the super- 
position principle in quantum mechanics should lead to 
the existence of a state that is thermalized with multiple 
temperatures. 
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